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ABSTRACT 

Airframe noise is a significant component of environmen- 
tal noise in the vicinity of airports. The noise associated with 
the leading-edge slat of typical transport aircraft is a prominent 
source of airframe noise. Previous work suggests that a slat-cove 
filler (SCF) may be an effective noise treatment. Hence, devel- 
opment and optimization of a practical slat- cove -filler structure 
is a priority. The objectives of this work are to optimize the de- 
sign of a functioning SCF which incorporates superelastic shape 
memory alloy (SMA) materials as flexures that permit the defor- 
mations involved in the configuration change. The goal of the 
optimization is to minimize the actuation force needed to retract 
the slat-SCF assembly while satisfying constraints on the maxi- 
mum SMA stress and on the SCF deflection under static aerody- 
namic pressure loads, while also satisfying the condition that the 
SCF self-deploy during slat extension. A finite element analysis 
model based on a physical bench-top model is created in Abaqus 
such that automated iterative analysis of the design could be per- 
formed. In order to achieve an optimized design, several design 
variables associated with the current SCF configuration are con- 
sidered, such as the thicknesses of SMA flexures and the dimen- 
sions of various components, SMA and conventional. Designs 
of experiment (DOE) are performed to investigate structural re- 
sponse to an aerodynamic pressure load and to slat retraction 
and deployment. DOE results are then used to inform the opti- 
mization process, which determines a design minimizing actua- 
tor forces while satisfying the required constraints. 


INTRODUCTION 

In transport-class aircraft, conventional high-lift systems 
(e.g., leading-edge slats and trailing-edge flaps) are used to aug- 
ment lift and improve stall characteristics at the low speeds re- 
quired for landing. These multi-element airfoil systems nest 
tightly in the cruise configuration to minimize drag. The de- 
ployed multi-element airfoil system, however, presents many ge- 
ometric discontinuities to the aerodynamic flow, such as edges, 
gaps and cavities. These geometric discontinuities are the cause 
for the unsteady aerodynamics that is the source for significant 
acoustic noise, termed airframe noise. The flow characteristics, 
noise production mechanisms and notional concepts for mitiga- 
tion of slat noise in particular have been studied extensively. The 
concept of the slat-cove filler (SCF) was introduced several years 
ago as a possible way to fill the cavity behind the deployed slat in 
order to reduce the unsteadiness in the flow and, thereby, reduce 
the radiated acoustic noise. Various idealized versions of the SCF 
concept have been considered by multiple research groups and 
the concept has been proven, both computationally and experi- 
mentally, to be effective at reducing slat noise [1-3]. 

This work considers to use of shape memory alloys (SMA) 
materials as flexures in a leading edge SCF that can be used to 
reduce the aeroacoustic noise produced by the wing of a typi- 
cal transport aircraft. Specifically, we perform analysis-driven 
design optimization of the SCF, which is currently in the devel- 
opmental stage. 

SMAs are attractive for use in morphing structures because 
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of the large recoverable deformation that can be achieved with- 
out incurring plasticity and their high energy density. Further, the 
use of an SMA component can reduce the number of parts and 
complexity of a system [4] . With regard to aerodynamic noise re- 
duction, the Boeing VGC [5] demonstrated how the shape mem- 
ory properties of SMAs can be used to achieve noise reduction. 
For this work, the superelastic properties of SMAs are exploited. 

Using the legacy method of smart structures design, obtain- 
ing an optimized solution required multiple iterations of design- 
ing, building, and testing physical designs. This resulted in a 
long and expensive cycle until an optimized solution is obtained. 
Due to the increasing capability of high-speed computational 
analysis, design problems can now be solved using optimization 
techniques in a more efficient way. Many of these techniques 
were discussed in a review of SMA design efforts [6]. Efforts 
have used genetic algorithms to explore the design space for ap- 
plications such as SMA spring actuation components [7], SMA 
wire-actuated rotors [8], and SMA structural damping mecha- 
nisms [9] . Another optimization technique is the gradient-based 
method which has been used for applications such as an active 
panel structure [10]. For this work, the goal of the optimization 
is to minimize the activation force needed to retract the slat-SCF 
assembly. 

This paper describes the concept of the SMA-based SCF and 
the process of its optimization. Specifically, Section 1 presents 
the physical form of the slat-SCF assembly. Section 2 addresses 
the computational framework and introduces the virtual model 
of the SCF. Section 3 discusses the goal of the optimization, the 
constraints, and the design parameters. Section 4 addresses the 
exploration of the design parameters and their influence on the 
model. Section 5 presents the results of the optimization. Sec- 
tion 6 summarizes the work and proposes concepts for future 
development. 

1 CONCEPT FOR SMA-BASED SLAT COVE FILLER 

As stated earlier, the goal of this work is to optimize the 
design of an SMA-based SCF. This proposed part fills the cavity 
behind the leading edge slat and guides the flow between the slat 
and the leading edge of the wing when the slat is deployed during 
low speed maneuvers [11]. 

Considering the amount of strain that will be needed for 
the SCF to be retracted (estimated to be about 2-5% [11]), typ- 
ical aerospace materials cannot be used because those materi- 
als would yield during the retraction of the SCF. In order to sat- 
isfy the deformation requirement, superelastic SMA components 
were considered. This type of material is assumed or chosen to 
be in the austenitic phase for all operating conditions and, when 
given a sufficient applied stress, will transform into the marten- 
sitic phase, allowing it to achieve large amounts of deformation 
without inducing plasticity. Once the load is removed, the SMA 
component will return to its original form; the deformation will 


be recovered as the material returns to the austenitic phase. 

In a related study, a physical 2-D prototype based on a typ- 
ical transport aircraft was developed [11]. The deployed and re- 
tracted configurations of the assembly can be seen in Fig. 1. The 
slat-SCF assembly consists of three parts; the leading edge slat, 
a hinge, and the SCF. The hinge is made of steel and is used to 
help facilitate the deployment and retraction of the SCF. The SCF 
can be made up of a forward SMA flexure, a steel mid-link, and 
an aft SMA flexure or can be formed from a monolithic strip of 
SMA. The slat-SCF assembly and the leading edge make up the 
overall leading edge assembly. When the slat-SCF assembly is 
in the process of retraction, the SCF will come into contact with 
the leading edge of the wing. This will cause an applied contact 
load between the SCF and the leading edge. Due to the proper- 
ties of superelastic SMAs, the leading edge will cause the SCF to 
deform into a configuration that can fit inside the confined space 
between the slat and leading edge during cruise. As the slat is de- 
ployed, the SMA flexures will lose contact (thus being unloaded) 
with the leading edge. If the configuration is correctly designed, 
this will result in the SMA components returning to their original 
form before retraction. 



Figure 1 . Picture of physical 2D model in its deployed configuration and 
its retracted configuration. Courtesy of [11]. 


2 ANALYSIS TOOLS and FEA MODEL 

With an introduction to the proposed SMA-based SCF pro- 
vided, the focus shifts to describing the computational frame- 
work that wis implemented and the analytical model of the phys- 
ical assembly. 

Computational Framework 

To perform the analysis of different configurations of the 
SCF and thus obtain an optimized solution, an efficient compu- 
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tational analysis framework is needed. The chosen framework 
uses custom constitutive model implementations in combination 
with a finite element suite. This combination is then automated 
by a simulation process management tool [6], which generates 
the desired finite element analysis (FEA) models and performs 
the analysis of those models. A post-processor is then used to 
obtain both local and global results of the analysis, which in turn 
can be used in a larger framework that assesses a design response 
and satisfaction of various design criteria and chooses new con- 
figurations to be analyzed. 

In order to analyze the SMA flexures, a constitutive model 
is needed that captures the necessary mechanical and thermal re- 
sponses. For this work, temperature is considered to be constant 
so the constitutive model must be able capture the superelastic 
effect of the SMA components. To do this, we have chosen 
the constitutive model developed by Auricchio and coworkers 
( [12], [13]) . The advantage of this model is that it has been pre- 
compiled into Abaqus, allowing team members to participate in 
the modeling who do not have the necessary compliers to run 
custom-coded user material subroutines [6]. 

This framework is flexible, allowing various FEA tools, con- 
stitutive models, or simulation process managers to be substi- 
tuted at any time. In this work, a combination of ModelCenter, 
Abaqus, and an Abaqus-native UMAT are used. ModelCenter 
is a design integration and optimization tool [14], which in this 
work directs the pre-processing, processing, and post-processing 
of Abaqus [15] through the use of custom scripts. 


FEA Model 

For this work, an FEA model based on the 2-D physical as- 
sembly is created and can be seen in Fig. 2. The hinge, slat, and 
leading edge are considered to be rigid. The SCF (a single ’’part” 
in Abaqus) is split into four sections; the SCF-hinge arm, the for- 
ward SMA flexure, the steel mid-link, and the aft SMA flexure. 
The partitions in the SCF are created using datum planes. This is 
done so that certain design parameters such as the lengths of the 
SMA flexures can be easily varied and the hinge arm can be ex- 
tended along the curve of the SCF without modifying the hinge 
part. The lengths of the various SCF sections are considered as 
the distance that a datum plane is offset from its reference plane. 

The mid-link is given material properties of steel (Young’s 
modulus of 2000ksi, Poisson’s ratio of 0.25, and density of 
0.000732 slinch/in 3 ) and a thickness of 0. 125in such that it would 
be rigid as compared to the SMA flexures. The SCF-hinge arm 
is also given material properties of steel. In order to implement 
the Auricchio model, mechanical constants (both material prop- 
erties and other parameters) must be specified. These constants 
are given in Table 1. The nominal values are based on typi- 
cal values for superelastic SMAs. The experimental values are 
based on actual hysteresis loops for this SMA material provided 
by Johnson-Matthey Inc. The density of the SMA material is 
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Figure 2. Analytical Model of the Leading Edge Assembly. 


taken as 0.000603 slinch/in 3 . If a monolithic SMA design is de- 
sired, the length of the mid-link can be reduced to a value that 
would essentially neglect the mid-link. 

In order to validate the experimental values, a hysteresis 
loop is generated in Abaqus using a two element SMA cube with 
the properties of Table 1 . The loop is then superimposed on the 
experimental hysteresis loop of the SMA material which can be 
seen in Fig. 3. Note that the test that obtained the experimental 
response did not start at zero strain so the computed hysteresis 
loop is shifted such that both responses started at the same strain. 
Note also that only the first curve is of interest. As seen in this 
figure, the hysteresis loop from Abaqus matches the experimen- 
tal response of the material quite well. 



Figure 3. Abaqus generated hysteresis loop superimposed on material 
response from an experiment. 
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Table 1. Nominal and experimental material properties used for the 
analysis of the SCF. Properties are associated with the Auricchio 
model [12], [13]. 


Constant 

Nominal 

Experimental 


(Elastic Properties) 

Ea 

7980 ksi 

7090 ksi 

Em 

6670 ksi 

5800 ksi 

* 

II 

0.33 

0.33 

(Phase Diagram Properties) 

0 Ms = a Ms 

59 ksi 

68.3 ksi 

a M f 

74.5 ksi 

75.1 ksi 

<3 As 

32.2 ksi 

28.5 ksi 

a A f 

21.5 ksi 

24.0 ksi 

% 

II 

1070 psi/ °C 

1450 psi/ °C** 


(Transformation Strain Properties) 


II 

5.6% 

4.4% 


(Other Constants) 

N a 

0 

0 

To 

300 

300 


*Nominal values; not experimentally obtained. 
** Assumed values; not experimentally obtained. 


There are two variations of the hinge and the SCF used in 
this study. The first variation is used in the analysis of the model 
response to aerodynamic loads. This variation has all three cut- 
ting planes placed horizontally with a horizontal reference plane 
at the forward end of the SCF. Each of the planes is offset from 
the plane below it in an additive scheme. For example, the ’’for- 
ward SMA” cutting plane is offset from the ”SCF-hinge arm” 
cutting plane while the latter is offset from the horizontal base 
plane. The first variation can be seen in Fig. 4. For this variation, 
approximately 840 linear shell elements (type S4R) are used and 
are evenly distributed throughout the SCF. To analyze the model 
response, static analysis (a Static, General step in Abaqus/CAE) 
is considered and the slat is stationary. 

Aerodynamic forces are considered while the SCF is fully 
deployed. It is assumed that these forces can be neglected when 
the slat-SCF assembly is in the process of being retracted or de- 
ployed. The aerodynamic forces are modeled as a distributed 
pressure load based on an X,Y,Z scatterplot (called a ’’point 



Figure 4. First Variation of the SCF and hinge. 

cloud” in Abaqus) that contained a dynamic pressure value at 
each of a number of points in 3-D space. This point cloud is cre- 
ated by linearly interpolating in the U1 (chordwise, X) and U2 
(airfoil thickness, Y) direction between data points that are com- 
putationally determined by NASA aerodynamicists. The pres- 
sure is taken to be uniform in the U3 (spanwise, Z) direction. 
Fig. 5 shows the pressure distribution (in psi) plotted along the 
SCF. Also for this first variation only, the nominal properties of 
Table 1 are used because the experimental values had not yet 
been determined at the time of this analysis. Given that this is 
an elastic analysis with small deflections, it is assumed that the 
results of the analysis would not vary significantly had the exper- 
imental values been used (see Ea and Em values in Table 1). 


Figure 5. Distributed aerodynamic pressure load plotted along the SCF. 
Pressure increments are in psi. 
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In order to consider configurations in which the hinge arm 
is smaller, the second variation of the SCF and hinge is used. 
For this variation, the arm of the hinge part is shortened while 
the forward end of the SCF is extended as seen in Figure 6. The 
”SCF-hinge arm” cutting plane is offset from a vertical reference 
plane that is positioned at forward end of the SCF. Also in this 
variation, the ’’forward SMA” cutting plane is now offset from 
the horizontal base plane instead of the ”SCF-hinge arm” cut- 
ting plane, and is given a minimum value such that ”SCF-hinge 
arm” cutting plane and the ’’forward SMA” cutting plane would 
not cross. The second variation of the hinge SCF is used in the 
retract-deploy analysis. 

In order to retract and deploy the slat-SCF assembly, a hinge 
connector (a feature in Abaqus which connects different bod- 
ies in an assembly) is implemented between a reference point 
(a specified point in Abaqus) and the hinge. The reference point 
is then given a rotational displacement in the U3 direction such 
that the slat-SCF assembly would pivot about the point. In order 
to fully retract and deploy the assembly, a displacement of 0.57 
radians is required. 

While the dynamic response of the slat/SCF system dur- 
ing retraction and deployment is not of primary interest (i.e., 
only the final states are assessed), the bi-stability and associated 
snap-through of the SCF necessitate the use of dynamic analy- 
sis methods. For purposes of efficiency and element flexibility, 
the Abaqus/Standard implicit dynamic solver is chosen. A single 
analysis of the retract-deploy case would take approximately 3-4 
hours using the previous variation. To make analysis more time 
efficient, the SCF is made to be one element thick with symmetry 
conditions in the U3 direction that are placed on both sides of the 
part, making this a ”2.5-D” analysis. To reduce the required run- 
time even further, the contact between SCF and the other parts 
is converted from infinitely stiff to a linear perturbation law (still 
very stiff). The use of the second variation of the hinge and SCF 
and the one element thick constraint results in a mesh with ap- 
proximately 140 linear shell elements (type S4). 

3 OPTIMIZATION PROBLEM 

Here we address the goals of the optimization, the con- 
straints that must be satisfied, and the design parameters that will 
be varied in order to obtain an optimized design. 

The goal of this optimization is to minimize the activation 
force needed to retract slat-SCF assembly. It is expected that a 
more compliant design will also generally correspond to a lighter 
design. This actuation force is measured in Abaqus as the ap- 
plied moment on the reference point from which the slat-SCF 
assembly pivots about. This is an important criterion because 
the smaller the actuation force, the less power that is needed to 
retract the entire slat, thus making this concept more viable for 
eventual implementation. A secondary goal is to minimize the 
contact pressure between the SCF and the leading edge. If the 
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Figure 6. Second Variation of the SCF and hinge. 


pressure becomes exceedingly large, the leading edge can be de- 
formed, possibly causing losses in desired aerodynamic perfor- 
mance of the leading edge. Any deformation of the leading edge 
would make the SCF not viable for use. The parameters that are 
varied consist of the length of the SCF-hinge arm (L/J, the length 
of the two SMA flexures^/y and Lf orwar( i ), and the thickness of 
the two SMA flexures (t a ft and tf orwar d). The desired shape of 
the deployed SCF is taken to be fixed (i.e. determined by aero- 
dynamicists). Therefore, one of the two SMA flexure lengths or 
the length of the mid-link can be directly influenced by changing 
the other two lengths. For this work, the aft SMA flexure is var- 
ied by varying the length of the forward SMA flexure (L f orwar d) 
and the length of the mid-link ( Lml )• 

The constraints under which this optimization are performed 
are that the SCF must self-deploy when the slat is rotated for- 
ward, the von Mises stress at any point in the SMA flexures does 
not exceed lOOksi during retraction, and that the SCF will be 
displaced no more than 0.1 in due to aerodynamic loads. If the 
SCF does not self-deploy, then it cannot modify the flow and is 
useless. Whether or not it self-deploys are determined by mea- 
suring the rotational displacement of the hinge. The stress that 
the SMA flexures experience must be constrained, because while 
the SMA components may be extremely flexible, it can still be 
plastically deformed due to excessive stress. If this occurs, then 
the desired effects of the SCF on the flow may not occur. The 
value of lOOksi was chosen because it is an approximate value of 
yield stress for this type of SMA material. The last constraint for 
this optimization is needed because when the slat-SCF assembly 
is deployed, the SCF will be subjected to aerodynamic loads. If 
the structure is displaced excessively, then its effect on the flow 
will be lessened. 
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4 INFLUENCE OF DESIGN PARAMETERS 

Before optimization could be conducted, it was necessary 
to perform a design of experiment (DOE) in order to understand 
how the design parameters effected the response of the structure 
to the various loads. It was decided that the aerodynamic load 
and the retract-deploy steps would be split into separate DOE 
studies so that model response to each case can be analyzed sep- 
arately. Due to the amount of information that is obtained from 
a DOE, only the most significant results are presented here. 

Preliminary Aerodynamic Load DOE 

The first DOE performed assessed how the SCF was affected 
by the aerodynamic load while the design parameters were var- 
ied. For this DOE, the simulation process control, FEA, and 
UMAT components of the analysis framework were used. With 
the Design of Experiment feature in ModelCenter, 243 combina- 
tions of L h , t afu t forward, L forward , and L ml were analyzed using 
a full factorial array. Out of the 243 runs, 60 combinations failed 
to converge. It was believed that these designs failed because 
they buckled which cannot be modeled via the approach used 
here. Had those runs converged, large values of displacement 
would have been observed. However, the remaining 183 suc- 
cessful runs allowed for valuable conclusions to be made about 
the response of the SCF to the aerodynamic load. The bounds for 
the parameters are shown in Table 2. 

Table 2. Design variables used in the exploration of the design space 
for the aerodynamic load DOE. L a f t ranges from approximately 1 .7in to 
10.1 in. 


Design Variable 

Range 

Lh 

0.2-0. 6 in 

L forward 

2-6 in 

Lml 

2-6 in 

taft 

0.02-0.06 in 

t forward 

0.015-0.06 in 


The main outputs from this analysis were the maximum Von 
Mises stress in the SMA flexures and the maximum magnitude of 
displacement anywhere on the SCF. These two parameters define 
the SCF response under the aerodynamic loading. 

A factor effects plot can be seen in Fig. 7 1 , where it was ob- 
served that the SCF thicknesses appear to have the most effect on 


^ote that the failure of many runs has resulted in an unbalanced DOE, lim- 
iting the utility of this plot in the strict sense. 


the maximum displacement while the length of the SCF-hinge 
arm has a negligible effect. More information about the rela- 
tion between the displacement and the thicknesses can be seen 
in Fig. 8, which shows a ModelCenter contour plot of the dis- 
placement in terms of the two most important variables. As seen 
in this figure, as either thickness increased, the maximum dis- 
placement decreased. However, if one of the thicknesses was 
small, the other must be thick enough to compensate for the loss 
in stiffness. Although it was not shown, another trend that was 
observed was that as the mid-link length increased, the maxi- 
mum displacement decreased. Clearly if more of the SCF was 
steel, the overall structure was stiffer, thus reducing the displace- 
ment. Of all the different combinations, no design corresponded 
above lOksi, well below the limit. 



Figure 7. Influence of the design parameters on the maximum displace- 
ment of the SCF. 


Retract-Deploy DOE 

The second DOE considered the results of many retract- 
deploy analyses. Like the previous DOE, all parts of the analysis 
framework were used. With the Design of Experiment feature in 
ModelCenter, 49 combinations of L h , t afu t forward, L forward , and 
Lml , using a Design Explorer Orthogonal Array, were analyzed. 
These 49 combinations were then combined with a smaller DOE 
set of 25 combinations. The bounds for the parameters for both 
DOE studies are shown in Table 3. 

The main model outputs of interest in this DOE were the 
rotational displacement of the hinge after deployment, the max- 
imum applied moment needed to retract the slat-SCF assembly, 
the maximum von Mises stress in the SMA flexures, and the max 
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0.06 in 



Figure 8. Contour plot showing the effect of the thicknesses on the max- 
imum displacement (L/ l =0.4in, Lf orwar d= 4in, Lml=4 in, L a f t ~ 5.9in). 

Table 3. Design variables used in the exploration of physical design 
space for the retract-deploy analysis. L a f t ranges from approximately 
2.3in to 12.2in. 


Design Variable 

Range 

Lh 

0.7-2 in 

L forward 

2-6 in 

Lml 

0.1-6 in 

Lft 

0.01-0.05 in 

t forward 

0.01-0.05 in 


contact pressure with the leading edge. Due to a computer fail- 
ure, the DOE with 49 combinations was disrupted so only 31 
designs were analyzed. Of these 31 combinations, 7 failed to 
converge. The DOE using 25 combinations also had 7 designs 
fail to converge. By assessing the output files associated with 
these runs, it was determined that analysis failure was due to a 
complex SCF deformation mode that is considered unacceptable. 
Therefore, failed runs are taken to indicate an infeasible design. 

The rotational displacement of the hinge after deployment 
was a measure of whether or not the SCF would redeploy au- 
tonomously. It was observed that this displacement was either 
less than 0.1 radians (SCF redeploys) or greater than 1 radian 
(SCF does not redeploy). It was also observed from the DOE 
results that when Lf orwar d and Lml were both greater than 4in, 
the SCF typically did not redeploy. Based on this, one can con- 


clude that the aft SMA flexure drived the redeployment. If it was 
not large enough, the aft flexure was not able to pull the SCF out 
of the slat cove and back into its original form. The region of 
the design space corresponding to a short aft flexure should be 
avoided during the optimization. 

Although not shown, another important result was that the 
length of the SCF-hinge arm had a significant effect on the max- 
imum moment. After a particular value of the hinge-SCF arm 
length (estimated to be between 1.35in and 1.65in), the mo- 
ment increased dramatically. This meant that large values for the 
length of the SCF-hinge arm should not be considered when op- 
timizing the design. A relation between the thicknesses and the 
maximum moment can be observed from Fig. 9, which shows a 
contour plot of the maximum moment. As seen in this figure, the 
lower values of maximum moment typically occured at a rela- 
tively small t f orwar d- Fower moments at smaller thicknesses was 
not surprising. Smaller thicknesses meant that the SCF was less 
stiff so there does not have to be as much applied moment to 
retract the assembly. However, based on the aerodynamic load 
DOE, excessively small thicknesses did not satisfy the displace- 
ment constraint. 



Figure 9. ModelCenter-generated contour plot showing the effect of 
the thicknesses on the maximum moment (L/*= 0.7in, Lf orwar d= 4in, 
Lml= 0.1 in, L a f t tt 10.2in. Generated from 42 total runs). 


One final relationship of interest was that of the relation be- 
tween the contact pressure and maximum moment. As seen in 
Fig. 10 the relation between the contact pressure and the max- 
imum moment was approximately linear. The more activation 
force that was applied to the slat-SCF assembly, the more load 
that was transferred from the SCF to the leading edge. As a re- 
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suit of this relation, optimizing the SCF for the maximum applied 
moment will also optimize the design for the contact pressure. 
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optimization, and designs can still satisfy the displacement con- 
straint with this upper bound. The lower bound on both t a f t and 
t forward was increased. Based on the DOE results of the pressure 
analysis, a thickness of 0.01 in was too thin, generally resulting 
in displacements greater than O.lin. 


Table 4. Design variable ranges used for optimization. L a f t ranges from 
approximately 4.3in to 12.2in. 


Design Variable 

Range 

Lh 

0.6-1. 5 in 

Lforward 

2-6 in 

Lml 

0.1-4 in 

taft 

0.02-0.05 in 

t forward 

0.02-0.05 in 


Figure 10. Relation between the maximum applied moment and the 
maximum contact pressure with the leading edge. 


5 OPTIMIZATION 

In order to set the bounds for the design parameters the re- 
sults from both DOE studies had to be considered. The chosen 
bounds can be seen in Table 4. Due to the large maximum mo- 
ment that results from having a large hinge arm, the upper bound 
of Lh was set to 1.5in. Changing this bound has negligible ef- 
fects on model response to the aerodynamic loading as Lh had a 
negligible effect on the displacement. Also, the lower bound of 
Lh was decreased so that lower values of the SCF-hinge arm can 
be considered. The new value for the lower bound was placed at 
the end of the arm for the hinge part. Since the length of the aft 
SMA flexure was an important factor in the redeployment of the 
SCF, the bounds of Lf orwar( i and Lml were considered at the same 
time. As stated earlier, when both lengths were above approxi- 
mately 4in, the SCF typically failed to redeploy. To avoid this 
region of the design space, the upper bound of Lml was set to 4 
in while the bounds of Lf orwarc i were not changed. The bounds of 
Lml were chosen to be adjusted because upon further inspection 
of the DOE results of the retract-deploy analysis, large values of 
Lml typically corresponded to large values of stress in the mid- 
link and a higher mass for the SCF. Also, larger values of Lml 
were more likely to result in the SCF not redeploying than larger 
values of Lf orwar d. Having a smaller upper bound for Lml will re- 
sult in designs with a larger maximum displacement. However, 
the constraint that the SCF redeploy is more important for this 


Based on the results of the retract-deploy DOE, it was ex- 
pected that the design that will minimize actuation moment will 
have the smallest thicknesses for the flexures. However, these 
thicknesses must not allow the SCF to fail the aerodynamic load- 
ing displacement constraint. Because of these conflicting trends, 
both the pressure analysis and the retract-deploy analysis were 
performed during the optimization. Both analyses were con- 
ducted using the same FEA model from the previous retract- 
deploy DOE. However the number of elements was increased 
to approximately 170 linear shell elements (type S4). This was 
done so that the number of elements along the edge of the SCF 
would roughly match the number of elements along the SCF edge 
that was used in the aerodynamic loads DOE. 

To perform the optimization, the Design Explorer tool in 
ModelCenter was chosen( [6], [14]). The optimization process 
required 283 FEA runs to converge to an optimized solution. The 
optimized values of the design input parameters and model out- 
puts are shown in Table 5. It is interesting to note that the op- 
timized solution occurred at the upper boundary for Lml • This 
means that a more optimized solution may exist if the steel mid- 
link was lengthened. However, as stated earlier, that region of the 
design space correlates to the SCF not fully redeploying. Future 
work will be needed to explore this region of the design space 
more intensively. 

6 SUMMARY AND CONCLUSIONS 

During low speed flight, typical transport aircraft deploy 
high-lift systems, such as leading edge slats, in order to the im- 
prove flight characteristics of the aircraft. However when de- 
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Table 5. Optimized design parameters and output 


Design Variable 

Value 

Nearest 

Bound 

SFC-hinge Length ( L h) 

0.647 in 

0.6 in 

Foward Flexure Length (Ef orwar d) 

3.28 in 

2 in 

Mid-Link Length (Lml) 

4.00 in 

4 in 

Aft Flexure Thickness (t a ft) 

0.04625 in 

0.05 in 

Foward Flexure Thickness (t forward) 

0.0238 in 

0.02 in 


Aft Flexure Length (L a yy)«7.02in 
Maximum Displacement^. 0991 in 


Rotational Displacements radians 
Maximum SMA von Mises Stress=70700psi 
Maximum Moment=913 lb-in 


ployed, these high-lift systems present geometric discontinuities 
to the flow. The implementation of a SCF is a proven means by 
which aeroacoustic noise reduction can be achieved. The over- 
all goal of this work was to optimize the design of a leading 
edge SMA-based SCF. The lengths of the hinge arm, forward 
SMA flexure, and mid-link and the thicknesses of the forward 
and aft SMA flexures were optimized in such a way that the acti- 
vation force, as measured by the applied moment, would be min- 
imized. The optimization was subject to the constraints that the 
von Mises stress in the SMA flexures would not exceed lOOksi, 
the SCF would redeploy autonomously, and that the SCF would 
not displace more than 0.1 in when subjected to given aerody- 
namic loads. A previously developed computational framework 
was used to perform, two DOE studies in order to investigate the 
model response to two different cases separately. The first DOE 
subjected the SCF to an aerodynamic load while it was deployed. 
The other DOE performed a retract-deploy analysis on the SCF. 
From the DOE studies, it was observed that optimizing the de- 
sign for either the activation force or the contact pressure with the 
leading edge would optimize the design for the other, and bounds 
for the optimization were chosen, avoiding regions of the design 
space where there would not be a solution. Using ModelCenter’s 
Design Explorer tool, both the aerodynamic load analysis and 
the retract-deploy analysis were performed in order to obtain an 
optimized solution. It was found that the optimial design occurs 
at the upper boundary of the mid-link length which implies that 
a more optimized solution may exist at larger values of Lml • 

The developments in this work have motivated future efforts 
such as comparing the analytical model to the physical model 


to in order to validate the results presented in this paper. With a 
validated model, further analysis such as the fluid- structure inter- 
action and in depth CFD can be conducted. Future work may im- 
plement other constitutive models so that a more accurate repre- 
sentation of the SMA material can be used. In addition, other de- 
signs for the SCF shall be considered so that aerodynamic noise 
can be reduced in the most optimum manner. 
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